Charm-sea Contribution to High-p T ip Production at the Fermilab Tevatron 
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The direct production of J/ip(4>') at large transverse momentum, pr ^> Mji$, at the Fermilab 
Tevatron is revisited. It is found that the sea-quark initiated processes dominate in the high-pr 
region within the framework of color-singlet model, which is not widely realized. We think this 
finding is enlightening for further investigation on the charmonium production mechanism. 

PACS numbers: 13.85. Ni, 12.38.Bx. 



Quarkonium production and decays have long been 
taken as an ideal means to investigate the nature of 
Quantum Chromo dynamics (Q CD) and other phenom- 
ena. Due to the approximately non-relativistic nature, 
the description of heavy quark and antiquark system 
stands as one of the simplest applications of QCD. The 
rich spectrum of its radial and orbital excitations pro- 
vides a suitable play ground for testing QCD based mod- 
els. The heavy, but not very heavy, quark mass enables 
one to get knowledge of both perturbative and nonper- 
turbative QCD via investigating quarkonium production 
and decays. The clean signals of quarkonium leptonic 
decays render the experiment detection with a high pre- 
cision, and therefore, quarkonium may play an unique 
role in the study of other phenomena as well, e.g. in 
detecting the parton distribution, the QGP signal, and 
even new physics. However, only with a theory which 
can precisely describe heavy quarkonium production and 
decays, may these advantages come true. 

During the past decade, intrigued by the discovery of 
J/ip(tp') surplus production at high pr at the Fermi- 
lab Tevatron [1,2], our understanding on the natures of 
quarkonium production and decays has experienced dra- 
matic changes. 

Conventionally, the so-called color-singlet model 
(CSM) was widely employed in the study of heavy 
quarkonium production and decays [3]. In CSM, it is 
assumed that the QQ pair produced in a high energy 
collision will bind to form a given quarkonium state only 
if the QQ pair is created in color-singlet state with the 
same quantum numbers as the produced bound states; as 
well, in the quarkonium decays the annihilating QQ pair 
will be in short distance and singlet with the same quan- 
tum numbers as its parent bound states. It is assumed in 
CSM that the production amplitudes can be factorized 
into short distance and long distance parts. The short 
distance sector is perturbative QCD applicable, while all 
the long distance nonperturbative effects are attributed 
to a single parameter, the wave function. That is, e.g., 
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The wavefunctions can be either determined phenomeno- 
logically through experiment measurement of quarko- 
nium leptonic decay rates, like 
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or calculated from potential models. 

CSM provides a prescription for calculations of not 
only the inclusive production rate of quarkonium states, 
but also their inclusive decay rates into light hadrons, 
leptons, and photons. Based on it many investigations 
had been carried out in past more than two decades, and 
at least qualitative description of quarkonium production 
and decays were achieved. Nevertheless, color-singlet fac- 
torization is only an ad hoc hypothesis. There are no 
general arguments to guarantee such a naive model may 
still work up to higher order radiative corrections. In fact 
people got know for quite long time it does not. Further- 
more, the attempt to incorporate relativistic corrections 
also met difficulties [4] . Although this kind of shortcom- 
ings of CSM had been known for many years, and some 
disagreements between theoretical predictions and exper- 
iment data existed, there was no big breakthrough in the- 
ory until the CDF group released [1] the data on large-pr 
J/tp production collected in the 1992-1993 run. The new 
data, which benefited from the advanced technology of 
vertex detector were free of the large background from 
B decays, put the CSM in an awkward situation, as the 
data differ very much from the leading order(LO) CSM 
predictions in both normalization and pt scaling. 

In 1993, Braaten and Yuan [5] noticed that at suf- 
ficiently high pt the dominant charmonium production 
mechanism is the production of a parton with large trans- 
verse momentum followed by its fragmenting into a char- 
monium state. With including the contribution by frag- 
mentation mechanism, the prompt J/ip data can be ex- 
plained within a small amount of error [6,7], where the 
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charm quark splitting into Xc states and then feeding 
down to J /ip> contributes overwhelmingly. Nevertheless, 
one can merely get a similar px asymptotic behavior 
for ip' production with the same scenario, and large dis- 
crepancy in normalization remained as well. This phe- 
nomenon was referred as the so-called " -0'-surplus" pro- 
duction or "-anomaly". 

In 1997, a measurement of direct J /ip production ex- 
posed [2], in which the higher excited states feeddown 
were stripped off. To one's surprise the new experiment 
result excesses the CSM prediction by a factor of ~ 30, 
the same as in ip' production. Nowadays, the former 
" ^'-surplus" problem turns to be the generic " ^-surplus" 
problem tcrminologically. 

A general factorization formalism [8] developed from 
the non-relativistic QCD(NRQCD) [9], which describes 
the inclusive heavy quarkonium production and decays, 
were established by Bodwin, Braaten and Lepage(BBL). 
NRQCD is formulated from first principles and the BBL 
approach allows relativistic and radiative corrections to 
be performed safely to any desired order. One of the 
striking advancements of the new development from CSM 
is that within the BBL framework the intermediate QQ 
state, which subsequently evolves into quarkonium states 
nonperturbatively, can be in both color-singlet and -octet 
configurations. At first order in v, the relative velocity of 
heavy quark, BBL and CSM are coincident in describing 
the S-wave quarkonium production. 

Based on the BBL formalism, Braaten and Flem- 
ing suggested to solve the ip' surplus production puz- 
zle via color-octet mechanism (COM) [10]. They pro- 
posed that the dominant ip' production at high px is 
through the fragmentation of a gluon into a cc pair in 
color-octet configuration, which will evolve into ip' non- 
perturbatively. Indeed they gave a well-fitted curve to 
the data and from which the non-perturbative matrix el- 
ement < Og ( 3 <Si) > was extracted with a magnitude 
being consistent with the estimation from NRQCD "ve- 
locity scaling rules" . After their pioneer work, hundreds 
of investigations have been performed in order to find 
either the signatures of color-octet states or any implica- 
tion of the new proposal to other phenomena [11]. 

The present situation is that on one hand the COM 
stands as the most plausible approach, up to now, in 
explaining the J/ip(ip') production "anomaly"; on the 
other hand, this scenario encounters some difficulties in 
confronting with other phenomena [11]. The most strik- 
ing crisis is the absence of high-py transversely polarized 
J/tp and ip' at the Tevatron in the first measurement from 
CDF [12]. According to NRQCD spin-symmetry, and the 
prescription that the dominant charmonium production 
mechanism at high pt is of a gluon splitting into a color- 
octet 3 Si charm quark pair, such polarized states should 
appear [13]. Therefore, to what degree the COM plays 
the role in quarkonium production is still an open ques- 



tion to my understanding. To find distinctive color-octet 
signatures and to eliminate the large errors remaining in 
different fits for corresponding matrix elements are cur- 
rently urgent tasks in this research realm, for both theory 
and experiment. 

In order to overcome the difficulties COM met, people 
tried to attribute large amount of high-pT events to in- 
trinsic transverse momentum of the interacting partons, 
suppose that the large uncertainties existed in the fc t - 
factorization are manageable and the k t would still man- 
ifest itself in not very small- a; [14]. To be noted that in 
the fc t -factorization formalism, the analyses suggest that 
the direct J/ip production is still dominated by color- 
octet contributions, but from and 3 -Fj 8 \ up to large 
transverse momenta of the order px < 20 GeV. 

Now that the difficulties for direct J/tp and its radial 
excitation ip' production at the Tevatron are the same 
within a small amount of error, as aforementioned, one 
may reasonably infer that the origins accounting for the 
large discrepancies between experimental data and the 
color-singlet description for both states would be the 
same. On this premise our investigation in this work 
will be restricted to J/tp for simplicity. The results and 
conclusions are applicable to ip>' . 

We notice that within the framework of collinear fac- 
torization, quarkonium production processes initiated by 
the sea-quark interactions have been paid less attention 
in previous investigations in color-singlet prescription. In 
Ref. [15] the sea initiated processes were considered, in- 
directly, for the large-pT J/ipi^') production, where the 
relative importance of valence and sea parton interact- 
ing processes was not distinguished. The leading or- 
der (LO) charm-sea interacting process was investigated 
in ref. [16], and found that it contributes negligiblely to 
the large-pT *P production as expected, since CDF data 
indicates that the da/dp T scaling favors fragmentation 
process behaving like l/p T which happens beyond the 
LO. Therefore, in order to estimate the sea quark con- 
tributions we need to consider the possible processes be- 
yond the LO. It is well-known that the source of the 
charm-sea distribution can be traced back to higher or- 
der gluon-gluon processes. And, because of the large log- 
arithmic term remaining in the NLO result of gluon-gluon 
to charm pair [17], the evolution effects of summing up 
all the large logs to any order tend to be important for 
the processes we are interested in here. Indeed, follow- 
ing evaluation shows that the naive NLO gluon-gluon 
to charm pair partonic process, which was considered in 
the calculations of [6,7], can not simply substitutes the 
charm-sea induced processes of our concern. 

On the other hand, although superficially the sea- 
quark interacting processes seems to be negligible, since 
the sea distribution probabilities are pretty small com- 
paring to those of valence quarks and gluons. Due to 
being in high energy and at high, but not very high, px 
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region, the s/i, the ratio of Mandelstam variables of the 
total center-of-mass energy squared in s and i-channels, 
kinematically suppresses the valence-quark and gluon ini- 
tiated processes relative to the sea interacting ones. And 
in the meantime high energy and large px enhance the 
sea quark densities inside the incident hadrons. To see 
this picture more clearly, let us have a close look at 
the fragmentation prescription for quarkonium produc- 
tion. Generally, quarkonium fragmentation production, 
A + B — > tjj + X, can be expressed as 

da(A + B -> H( PT ) + X) = V / dx a f a/A (x a ) 

a,b,c J ° 

x / dx b f b/B (x b ) / dzda(a + b — > c(p T /z) + X) 
Jo Jo 

x D c ^ H (z, h) , (3) 

where c is the fragmenting parton, cither a gluon or a 
(anit)charm quark, and the sum runs over all possible 
partons. D(z,^i) is the fragmentation function and z is 
the momentum fraction of the the fragmenting parton 
carried by quarkonium state. The evolution of the frag- 
mentation function D c ^h(z, m) with scale \i in Eq.(3) 
is accomplished by the utilization of Alterelli-Parisi(AP) 
equations 

d f 1 dy 

°M • Jz V 

where the are the splitting functions of a parton j 
into a parton i. 

To show the importance of the sea quark interacting 
processes in the fragmentation approach (3), we do a 
simple comparison, for example, of the hard-scattering 
processes 

g + g^C + C (5) 

with 

g + C(C) - g + C(C) , (6) 

where the C and C stand for charm and anticharm 
quarks, which are produced slightly off-shell and in high 
energy with large px- In LO and massless limit, the dif- 
ferential cross sections for processes (5) and (6) are: 

and 

respectively. Since we are interested in large-p^ J/ipiip') 
hadroproduction at high energy, obviously in a certain 



scope of phase space, the process (5) is suppressed with 
respect to process (6) by the factor of s/t. However, this 
is just a schematic argument. To be more strict, we need 
to convolute the hard scattering cross section with parton 
distribution and fragmentation functions. Without losing 
qualitative correctness, for simplicity we do a comparison 
for the subsets g + g — > C + C and g + C(C) -> g + C(C) 
induced processes in pp collision by convoluting the hard 
part with only the parton densities, or in other words, 
by integrating out the common fragmentation probabil- 
ity J Dq . Direct numerical calculation shows that the 
cross section induced by hard process (6) overtakes what 
induced by (5) by a factor of 7 at px — 15 GeV. 
For completeness, we consider processes 

q l (q i ) + C(C)^q l (q i )+ C(C) , (9) 

in our numerical calculation as well. Here, the 
qi(q~i) represents partons of both valence and sea 
quarks(antiquarks) of the colliding nucleons. Practical 
exercise, similar as performed in preceding paragraph, 
shows that this kind of hard interaction processes also 
contributes more to J/ip large-pr production than via 
subprocess g + g — > C + C. And, to be noted that the 
latter was taken to be the dominant process in many of 
previous analyses within CSM. 

In the numerical calculation of the differential cross 
section, we need to choose a set of parton distributions. 
In this work we take CTEQ5M [18] parameterization as 
our input. We have also tried another set of parton dis- 
tributions, the MRST99 [19], and found that different 
parton distribution functions give similar results within 
tens of percent. That means the conclusions given in 
this paper will not be spoiled by taking a different set of 
parton distributions for convolution in Eq.(4). 
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FIG. 1. Fragmentation Contributions to the differential 
cross section for direct J/tp production at the Fermilab Teva- 
tron, compared with the CDF experiment data read from [2] . 
The upper curve corresponds to fiR — hf = /Ltfrag = Pt/2, 
while the lower curve to fiR = /Uf = AWrag = 2pr- 
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In figure 1, various fragmentation contributions to di- 
rect J/ip production are shown, and the sum of them 
is confronted to the CDF data. Since we merely want 
to sketch the importance of sea-quark contributions to 
the large-pT charmonium production, the lowest order 
hard scattering cross sections and fragmentation func- 
tion are employed. We also neglect the effect from non- 
diagonal splitting function, P cg , in Eq. (4), which, as 
pointed out in Rcf. [6,7], may give a factor of as large as 
1.5 to the charm fragmentation process. We estimate 
the uncertainties of higher order corrections by vary- 
ing the scale, given that a lower value of scale will ac- 
count for, in a certain degree, the contributions from 
higher order corrections and non-diagonal splitting func- 
tions. The upper and lower lines in the figure are ob- 
tained by varying the scales of factorization, renormal- 
ization and fragmentation. The upper curve corresponds 
to pn = jip = Mfrag = Pt/2, while the lower curve to 
A*-R = (J'F = Mfrag = %Pt- In drawing the diagram we 
use the fragmentation function given in Ref. [5] and val- 
ues quoted thereof (Rq = 0.8GeV 3 , a s = 0.26, m c = 1.5 
GeV) . The symmetry of a sea quark and its antiquark in 
hard scattering and fragmentation precesses is invoked. 
In addition, we perform our calculation in the pr > 10 
GeV region, where contributions from lowest order par- 
ton fusion processes can be safely neglected, 

From figure 1 we see that the discrepancy between di- 
rect J ftp production data and the CSM prediction(the 
sum of different processes) is less than an order in optimal 
case (the upper solid line) , rather than the common belief 
of 30 or more. The finding that the charm-sea initiated 
processes contribute dominantly to high pr charmonium 
production, as shown in figure 1, within CSM looks sur- 
prise, whereas, it is not really an unthinkable thing. Simi- 
lar cases exist in some other quarkonium high energy pro- 
duction processes as well. For example, in photon-photon 
collision, quarkonium production via resolved processes 
is not always minor to the direct one. In addition, it 
should be noted that although the charm sea originates 
from the high order valence quark and gluon interactions, 
the naive NLO QCD result for J/ip(ip') hadroproduction 
can not simply substitutes the result from sea-qaurk ini- 
tiated processes. In the latter case we are considering 
the production via fragmentation mechanism, which is 
already beyond LO in a s . 

To conclude, in this work we study the relative impor- 
tance of the sea quark initiated processes with respect 
to the valence quark and gluon initiated ones for large- 
Pt J/tpi^') production within the CSM. It is found that 
the former may contribute more than the latter by a fac- 
tor of six. We notice that to many people within the 
community the gluon initiated processes are still taken 
to be the dominant ones for the high-pr J/ip(ip r ) pro- 
duction from CSM calculations. We hope this work may 
elucidate it somehow. With including the new produc- 
tion scheme, the total cross section from CSM prediction 



fall off the experiment data by less than an order in ex- 
treme situation. Therefore, to explain the CDF data, 
COM is still necessary and essential. Nevertheless, the 
increase of the contribution from CSM means a shrinkage 
of the contribution from COM, which might be as large 
as twenty percent. Furthermore, the result in this work 
gives us a strong hint that the J/tp(ip')-smplus produc- 
tion at the Fermilab Tevatron might still be explained 
within CSM after including the NLO calculation, like 
the charmonium photoproduction at HERA [20] [21]. At 
least the color-octet contribution will be much less than 
previously thought. Finally, we would like to point out 
that since the (anti)charm quark is taken to be massless 
in our calculation, some uncertainties will be induced by 
this measure to order of 0{Mj/^/pt)- For detailed esti- 
mations of the uncertainties by taking zero mass scheme 
see Refs. [22] [23]. In all, the error would be within a 
factor of two in our analysis, which will not change the 
conclusions of this work. 

This work was supported by the Grant-in- Aid of JSPS 
committee. The author would like to thank K. Hagiwara 
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